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INTRODUCTION  AND  SUMMARY 


In  designing  guidance  or  detection  systems,  an  important  aid 
towards  optimal  or  quasi-optimal  design  is  a  knowledge  of  the  proper¬ 
ties  of  the  background  in  which  the  system  must  work.  The  present 
paper  is  concerned  with  the  effect  of  background  gradients  on  infrai  ed 
systems  that  use  a  rotating  reticle  in  the  focal  plane  of  a  scanning  sys¬ 
tem  AS  an  aid  in  discriminating  against  background  signals  in  favor  of 
desired  or  target  signals. 

The  background  radiation  is  regarded  here  as  a  two-dimensional 
random  function  of  space  or  angular  variables.  (The  idea  of  treating 
the  background  in  this  manner  is  not  original  with  the  writer  but  has 
appeared  in  the  works  of  R.  C,  Jones  and  perhaps  others.  The  writer 
knows  of  much  important  work,  by  Dr.  Jones  in  this  field. ) 

The  scanning  and  rotating  action  may  be  said  to  convert  the 
background  or  input  random  process  or  "noise"  into  a  new  noise  pro¬ 
cess  at  the  output  side  of  the  reticle.  If  this  output  noise  can  be  de¬ 
scribed  in  terms  of  an  output  correlation  function  or  pow'er  spectrum, 
it  becomes  possible,  in  a  limited  sense,  to  calculate  the  effect  of  pro- 
cessixig  by  photoconductive  cells,  amplifiers  and  the  like.  Thus,  the 
system  analysis  stage  of  design  may  be  facilitated  and  possibly  a  little 
prioress  made  toward  more  sophisticated  solutions  of  the  design  syn¬ 
thesis  problem. 

In  this  paper,  expressions  are  derived  for  the  output  autocorre¬ 
lation  quantity^y(t)  y(t  where  y(t)  and  y(t  +  t)  are  the  reticle  out¬ 

puts  at  time  t  and  t  +  t  reifpectively .  and  ^  prefers  to  an  ensemble 
average.  If  the  background  random  function  may  be  regarded  as  a 
Stationary  random  function,  this  autocorrelation  function  may  be  con¬ 
verted  by  JTourier  transformation  to  an  output  power  spectrum. 

The  analysis,  which  is  thought  to  contain  some  new  results,  has 
tihe  advamtage  of  considerable  generality  and,  at  the  same  time,  has 
SWffifisTit  flexibility  to  allow  computation  in  simpler  cases  The  work 

also  to  clarify  the  roles  of  the  pertinent  variables  in  the  over-all 
If  fullest  generality  is  maintained,  the  expressions  obtained 


can  be  converted  into  numbers  only  with  the  aid  of  very  extensive  com¬ 
putational  aid.  In  certain  cases,  however,  the  expressions  may  be 
reduced  to  relatively  illuminating  forms,  and  the  computation  required 
may  be  reduced  to  a  moderate  amount. 

The  form  of  some  of  the  results  obtained  serves  to  point  up  the 
type  of  background  measurements  that  would  be  necessary  for  more 
accurate  estimates  of  the  outputs.  The  autocorrelation  function  of  the 
background  enters,  in  a  fundamental  way,  into  the  expressions;  it  may 
therefore  be  advisable  to  design  the  measurement  program  to  measure 
directly  this  correlation  function.  For  some  background  "bucking'*  or 
"suppression"  schemes,  success  in  cancelling  out  background  depends 
ultimately  upon  the  existence  of  relations  or,  somewhat  less  generally, 
of  cross -correlations  between  the  random  processes  in  the  various 
"colors.  "  This  would  suggest  the  investigation  of  a  program  of  measur¬ 
ing  these  cross -correlations. 

A  separate  study  (Appendix  B)  is  made  of  the  problem  of  studying 
background  by  a  technique  of  circular  scanning;  that  is,  scanning  the 
background  with  a  small  field  of  view  that  moves  repeatedly  around  a 
circle,  and  relating  the  Fourier  coefficients  of  the  resulting  periodic 
function  to  the  power  spectral  or  correlation  properties  of  the  back¬ 
ground,  An  analysis  is  made  of  the  processes  required  to  effect  this 
conversion  of  the  measurement  data.  The  limitations  of  this  procedure 
are  such  as  to  emphasize  the  desirability  of  having  direct  measurements 
of  the  correlation  function. 


2 


BASIC  FORMULATION 


Consider  an  optical  system  with  a  rotating  reticle  of  radius  a,  in 
the  focal  plane,  scanning  at  time  t  a  portion  of  the  background  in  a  cer¬ 
tain  region  of  the  spectrum  e.g.,  the  neighborhood  of  Z  microns.  Use  a 
set  of  polar  coordinates  (r,9)  {Figure  1)  to  describe  position  on  the  reticle. 
Since  that  portion  of  the  background  in  the  instantaneous  field  of  view'  is 
imaged  in  the  reticle  plane,  we  may  also  consider  r,  9  to  be  the  angular 
p<^ar  coordinatea  of  a  point  in  the  background.  Let  the  transmission 
function  of  the  reticle  be  R(r,9).  and  the  radiance  function  of  the  back¬ 
ground  to  be  B{r,9).  Then  the  expression 


y{t) 


-L 


a  Zn 

dr  d9rB(r.9)  R{r,e)  , 

0  Jo 


(1) 


may  he  regarded  as  proportional  to  the  output  of  the  reticle  at  time  t. 

The  factors  of  proportionality  omitted  involve  the  size  of  the  collecting 
system  and  the  transmission  factor  of  the  optical  system.  Henceforth 
we  shall  refer  to  y{t)  simply  as  the  "output. " 

Let  the  scanning  system  move  with  velocity  v  in  the  direction  of 
the  reference  angle  9  =  0  (Figure  1)  so  that  at  time  t  +  t  the  center  of  the 
Bsdd  of  view  has  moved  an  angular  distance  vt  in  the  direction  of  9  =  0. 
Let  the  coordinate  system  translate  with  this  scanning  motion  and  denote 
tite  coordinate  system  for  the  translated  reticle  by  (r',9')*  The  output  at 
time  t  4  T  is 


y<t  4  t) 


2it 


de'  r’  B(r',0')  R(r',9'  -  ut)  , 


(2) 


T - 

F<nr  this  to  be  strictly  true  B(r,9}  should  not  be  the  true  radiance 
ieactioiit  but  the  radiance  function  as  modified  by  the  aberration  of  the 
■yatam.  The  analysis  applies,  strictly  speaking,  to  this  modified  defi- 
aittop  of  B(r,0).  If  the  aberration  is  small  the  two  quantities  are  almost 

Mme. 

ibnotlior  modification  is  necessary  if  the  image  of  the  field  of  view  does 
Bat  coincide  with  the  reticle  circle,  or  if  there  is  another  significant 
^jasmiocion  function  of  (r,9)  in  addition  to  that  supplied  by  the  reticle. 
Vaa  aach  a  case,  the  work  through  Equation  (4)  would  be  unchanged  ex- 
aiit  ttet  B(r.«)  would  be  multiplied  by  an  "aperture  function"  say  A(r,9). 
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ROTATING  AND  TRANSLATING  SYSTEM. 

THE  SHADED  SECTION  INDICATES  A  RUDIMENTARY  RETICLE  FUNCTION. 


since  the  reticle,  which  is  rotating  with  angular  %'^elocity  w,  has  rotated 
through  the  angle  wt  while  translating  through  the  distance  vt. 

Let  us  form  the  product 


y(t)yCt  +  t) 


I  I 

-^0  Jo 


dr ' 


r 

-^0 


Zir 


de 


f 

■^0 


,-2ir 


d0’  rr*R{r,9)R(r',8' 


WT)B{r ,  9)B(r 8' 

(3) 


Now  the  ensemble  average  of  y(t)  y(t  +  t),  that  is  the  value  of 
yit)  y(t  +  t)  averaged  over  various  samples  of  the  background,  is 


^<t)y(t  +  J 


a 


rZrt 


rlf* 


dr 


d0 


0 


-0 


-^0 


2it 


d0'  rr'  RCr.BlRCr'.e-  >  «T)<B{r ,  0)  B(r ' ,  0 

(4) 


where  the  brackets mean  ensemble  average.  The  quantity  <^y(t)y(t  +  t|> 
is  the  autocorrelation  function  of  the  reticle  output  y{t).  If  the  background 
random  function  B<r,0)  may  be  regarded  as  a  "stationary"  random  function 
that  ia,  roughly  speaking,  if  its  average  properties  do  not  vary  with^  (r,0) 
over  a  region  interest,  then  we  can  define  a  useful  "power  spectral 
denaity”  or  "power  spectrum"  which  is  the  Fourier  transform  of 
V^fhrtt  +  .  The  effect  of  subsequent  processing  upon  the  reticle  output 

may  then  be  eaaUy  computed  in  a  limited  statistical  sense. 


The  central  goal  of  our  analysis  is  then  to  evaluate  <;^y(t)y(t  +  t)^ 
from  Equatiem  (4).  The  background  quantity  that  we  need  to  compute  this 
la  ^B(r,d)B{r',0*^  which  is  the  correlation  function  of  the  background. 
How,  in  particular  cases,  even  if  B|r.0)  may  be  regarded  as  a  stationary 
function,  it  may  happen  that  <^B(r.  0)  B{r ',  6')^  must  be  regarded  as  a 
function  of  the  vector  distance  between  the  points  (r.  0)  and  (r',  0').  If  this 


'€»  a  more  detailed  mathematical  description  of  a  stationary  random 
>M«aa  the  reader  is  referred  to  the  extensive  mathematical  literature 
ma  aubject  of  random  functions. 
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is  true  then  the  evaluation  of  Equation  (4)  will  probably  be  extremely 
difficult  (but  still  possible  with  the  aid  of  the  high-speed  digital  com¬ 
puters).  This  case  will  ..ot  be  considered  further  here. 


ADDITIONAL  ASSUMPTIONS 


in  order  to  make  the  general  result  of  Equation  (4)  more  amen¬ 
able  to  further  analysis,  we  shall  make  some  simplifying  assumptions. 

An  assumption  that  appears  to  be  useful,  is  that  the  autocorrelation  func¬ 
tion  <^{r,  6}  B(r 6*)^is  a  function yfd)  only  of  the  scalar  distance  d 
between  (r,  0)  and  (r’,8');  in  this  case  we  speak  of  the  background  random 
function  as  being  isotropic.  With  the  aid  of  Figure  1  it  is  seen  that  the 
distance  d  is  given  by 

/  2  2 
d  =  (r'cosS  +  VT  -  r  cos  0)  +  (r'sin©'  -  rsinB) 

=  r'^  +  -  2rr'cos{6'  -  9)  +  2vT(r'cos8'  -  r  cos  9)  +  (vt)^  .  (5) 


With  this  assumption,  Equation  (4)  takes  the  form 

<y(t)y(t  +  t)>  =  /  I  }  drdr'ded0'rr'R(r,e)R{r'.e'  -  WT)-)(r<d)  .  {4a) 


Another  simplifying,  but  useful,  assumption  is  that  R(r,9)  is  a 
function  only  of  9.  There  are  cases  when  reasonably  simple  results  can 
be  obtained  for  more  general  functions  R(r,8),  but  our  illustrative  ex¬ 
amples  will  be  limited  to  the  simpler  case.  Another  special  assumption 
of  useful  scope  is  that  the  scanning  speed  v  is  small  compared  with  the 
speed  of  most  points  of  the  rotating  reticle,  so  that  the  error  introduced 
in  setting  v  =  0  in  Equation  (5)  is  not  serious.  The  important  simplifi¬ 
cation  resulting  from  letting  v  =  0  is  that  d  becomes  a  function  of  0'  -  9. 


Incidentally,  for  the  case  in  which  the  field  of  view  is  rectangular  and 
the  reticle  motion  consists  of  parallel  bars  sweeping  across  the  field  of 
view,  the  assumption  v  =  0  is  not  necessary;  i.e,,  results  similar  to 
those  of  the  subsequent  example  may  be  obtained  even  if  v  ^  0. 
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A  CASE  OF  INTEREST 


With  these  assumptions  and  the  notation  Q(t)  . 

jEk^uation  (4)  becomes 


rZft 

rZv 

/o  “‘J 

1  de*  rr'  R(9)  R(0'  -  wt)! 

0  f 

where 


=  ♦  r*^  -  Ztt*  cO8(0'  -  9) 


1/2 


(6) 


(7) 


The  expression  (6)  has  some  interesting  properties.  First,  let 
US  rewrite  it  as 


■s: 


*2w  rZn 

Qtr)  -\  d«  /  de*  G(e'  -  6)  R{6)  K{^'  -  wt)  , 

Jo 


(8) 


where 


G<«’  -  9)  drj^  dr*  rr 


^  +  r^  -  2rr'  cosCe*  -  6)j 


(9) 


It  will  be  convenient  to  write  u  for  0'  -  6,  giving 
G(u)  s  f  dr  f  dr*  rr'H  fr'*"  +  r^  -  2r  r*  cos  ul 

Jo  Jo  )  I  J 


1/2 


(9a) 


An  important  result  relative  to  Equation  (8)  may  now  be  stated: 
Q(r)  may  be  written  as 


Q(t) 


rZJt 

=  I  du  G(u  +  wt)  P(u)  , 

Jo 


(10) 


where 


P«tt) 


-'0 


d9  R(e)  R{e  4  u)  . 


(11) 


The  method  by  which  this  result  was  obtained  is  given  in  detail  in  Appen¬ 
dix  A.  The  "intuitive"  significance  of  Equation  (iO)  and  (11)  is  as  follows*. 
The  form  of  Equation  (11)  indicates  that  P(u)  is  a  type  of  correlation 
function  to  be  called  hereafter  the  "reticle  correlation  function."  Q(t)  is 
the  correlation  function  of  the  oiftput  as  previously  discussed.  Equation 
(10)  expresses  Q(t)  as  a  convolution  of  G{u)  and  P(u}.  We  may  therefore 
regard  G(u)  as  another  correlation  function.  It  is  the  correlation  function 
of  the  background  as  "seen"  by  the  rotating  reticle;  it  may  be  called  the 

4 

"noise  correlation  function.  " 

Since  Q(t}  is  an  autocorrelation  function,  its  Fourier  transform 
is  a  "power  spectrum.  "  Taking  the  transform  of  Equation  (10)  we  have  an 
example  of  the  well-known  result  that  the  Fourier  transform  of  a  corre¬ 
lation  integral  is  the  product  of  the  Fourier  transform  of  the  convolved 
members. 

Of  course,  since  P(u)  and  G(u)  are  both  periodic  in  u,  Q(t)  is 
also  periodic  with  the  time  period  Zn/u.  Thus  the  power  spectrum  is  a 
line  spectrum  at  frequencies  which  are  multiples  of  w/Zf(. 

The  Fourier  transform  of  Equation  (10)  is  found  without  difficulty 

to  be 


where 


W  (f)  =  1  W*  (f/w)  W  (f/w)  W*  (f/w)  ; 

W 


W^..)  =  I 

w^,.,  =  / 


V  -Zirizu  J 
G(u)e  du  , 


-^TTiZU  . 

R(u)e  du  , 


If  V  had  not  been  set  equal  to  zero,  the  noise  correlation  function  would 
not  exist  in  this  form,  namely  as  a  function  of  u  =  0*  -  0,  For  the  rectang¬ 
ular  geometry  noted  in  footnote  3,  the  noise  correlation  function  exists  as 
a  function  of  x*  x  +  vt  where  x*  and  x  are  coordinates  which  are  rectang 
uiar  analogues  of  8^  and  9,  and  vt  has  the  same  meaning  as  above. 
Relations  similar  to  Equation  (10)  and  (11)  also  hold  for  this  rectangular 
case. 
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and 


Wq(£| 


Q(t)« 


-  Zirifr 


dT  . 


(15) 


The  eyxnhoi  denotes  the  complex  conjugate. 

Since  G(u)  is  an  even  function  of  u.  W^^z)  =  Wq(z),  and  we  may 
as  well  write 


(16) 


Since  G(u),  R{u),  and  Q(t)  are  periodic  functions,  the  definitions 
of  Equations  ( 13),  (14),  and  { 15 )  need  to  be  properly  interpreted.  Using 
a  delta  function  notation,  it  is  easy  to  show  that 


where 


where 


and 


where 


Wo(z) 

«g'*' 

W^(.» 

«R<*> 

V*' 

«Q<0 


oo 

23  5(z-n/2ir)  Uq(z)  . 

n=-co 

l/2ir  ^  G(u)e  du  , 

Jo 


00 

6(z  -  n/2x)  Uj^(z) 


,  ,  -2-irizu  j 

R(u)e  du 


oo 

^  6(f  -  nw/Sx)  U  (f) 

n=  -  oo 

/>  -Ztrifr 

w/2ir  Q(T)e 

Jo 


{13a) 


(13b) 


(14a) 


{14b) 


{15a) 


dr  . 


(15b) 


Thus  is  ^  delta  function  or  line  spectrum.  It  is  non-asero 

only  for  integral  multiples  of  at/Zir  .  The  *4ine-*  character  of  the  spectrum 
is  due  to  the  approximation  v=0.  If  v  is  slightly  different  from  sere, 
the  lines  are  changed  to  narrow  bands. 

Equation  (16)  or  the  equivalent  equation  (10)  exhibits  the  power 
spectrum  of  the  output  as  a  product  of  a  -^reticle  spectrum’-  and  a  ’’noise 
spectrum’’  where  the  noise  spectrtim  is  the  Fourier  transform  of  G(u). 

We  see  then  that  once  G(u)  or  its  transform  is  computed,  the  output  spec¬ 
trum  may  be  found  either  in  terms  of  the  convolution  integral  (10)  or  the 
spectral  product  (16).  The  "reticle  spectrum"  "reticle 

correlation  function"  P{u),  is  easy  to  obtain  in  terms  of  the  reticle  geo¬ 
metry,  and  the  definitions  of  Equations  (13)  to  {15b).  Thus,  once  G(u) 
is  found,  the  output  for  various  reticle  functions  is  simple  to  compute. 

The  problem  of  obtaining  G(u),  or  its  transform  Wq{z)  by  measurements 
or  calculations,  is  therefore,  of  key  importance.  This  problem  we  now 
examine. 

DETERMINATION  OF  NOISE  CORRELATION 
FUNCTION  OR  POWER  SPECTRUM 

From  Equation  (16)  one  sees  that  by  measurements  of  ^^(f),  the 
power  transform  of  Q(t),  and  measurements  or  calculations  of  W^(z),  one 
can  be  division  obtain  Wq(z),  the  transform  of  G(u}.  This,  in  principle 
is  a  valid  method  of  obtaining  Wq(2),  and  indeed  may  turn  out  to  be  the 
most  practical  method.  The  reticle  function  R(9)  may  be  chosen  for  con¬ 
venience.  Reticle  functions  will  generally  be  such  that  Wj^(z)  is  zero  for 
certain  values  of  z,  so  several  reticle  functions  may  have  to  be  used  to 
get  enough  information-  One  may  also  consider  using  a  "delta  function" 
reticle  (single  narrow  transnrutting  slit)  to  pass  "all"  values  of  z.  The 
results  obtained  for  Wq(z)  will  be  valid  only  for  the  particular  value  of  a 
used  in  the  measurements,  as  reference  to  Equation  (9a)  shows.  Thus 
one  would  havf*  to  repeat  the  process  for  each  value  of  a  considered. 

Another,  somewhat  more  fundamental  approach  is  possible  which 
we  now  wish  to  study. 
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JLet  us  consider  the  noise  correlation  fiuiction  G(u).  as  defined  in 
liquation  (9a}, 


G(u)  =  C  f"  dr  dr’  rr'/' 

Jo  Jo 


Z  2 

r  +  r*  -  2  rr’  cos 


“] 


i/2 


(9a) 


it. 


From  the  form  of  this  integral  we  see  that  G(u)  is  determined  by 
tile  values  of  ^{d)  for  the  range  of  d  between  0  and  2a.  (If  v  is  not 
aero,  reference  to  Equation  (4a)  or  Figure  t  shows  that  values  of  ^  tor 
a  aomewhat  larger  upper  value  of  d  are  needed  for  the  computation  of 
Q(t}.  This  suggests  that  direct  measurements  for  such  ranges 

{ratihier  than  measurements  such  as  have  been  made  up  to  the  present,  of 
tiha  power  spectrum)  might  be  undertaken.  For  the  purpose  of  obtaining, 
tltrough  expressions  such  as  have  been  derived,  insight  into  the  output 
neiae  of  the  reticle,  the  correlation  fiuiction  is  the  quantity  most  directly 
needed.  In  principle  one  can  Fourier -transform  spectral  measurements 
lato  correlation  information,  but  in  practice  large  errors  may  be  intro¬ 
duced  in  attempting  to  transform  the  necessary  incomplete  measurements. 

To  proceed  further  on  the  theoretical  side,  we  may  try  to  find 
ll^lauaible  forma  for  ^(d),  which  may  be  in  agreement  with  the  measure - 
 s  vit*  ^gaxita  ao  far,  and  see  how  these  influence  G(u)  aund  its  transform.  As 

i^geMhstioned  above,  the  limited  amount  of  available  experimental  material 
Iwriiig  on  this  question  has  been  expressed  in  terms  not  of  ^ ,  but  of  its 
 transform.  Thus  our  comparisons  at  present  between  theoreti- 
4Biikt  u^rsssions  and  experiment  can  be  made  most  conveniently  in  terms 
power  spectra*  The  two-dimensional  power  spectrum 
tbe  two^^dimensional  correlation  function  y)  are  Fourier  trans- 

ililmis  of  each  other.  That  is* 


xr  xr 


7Ti(k  X  +  k  y) 
X  y’ 


dx  dy  , 


(17) 


Xoo  ^OO 

j  * 

oo  “'-ao 


*^y*  2iTi(k^x  +  k^y)  j  dk^  dk^  . 


(I8j 


The  writings  and  measurement#  of  R.  C*  Jones  and  others,  including 
mcasureinents  of  sky  background  at  The  Ramo- Wooldridge  Corporation, 
suggest  the  form 


W^(k  ,  k  ) 

B  X  y' 


M 


(i9) 


where  M  is  a  constant,  at  least  for  f<.irly  high  wave  numbers  (number  of 
waves  per  radian).  More  exactly,  the  measurements  so  far  are  not  incon¬ 
sistent  w'ith  this  fsorm.  However,  tliis  form  is  not  acceptable  for  low  wave 
numbers,  because  the  integral 


dk 


X 


(20) 


that  is,  the  integral  of  the  spectrum  over  all  wave  numbers,  diverges- 

The  infinite  value  is  due  to  the  pole  at  k  =  k  =0. 

X  y 

A  way  to  adjust  Equation  (19)  so  as  to  eliminate  this  trouble  and 
still  not  be  inconsistent  with  the  experimental  data,  has  already  been 
suggested  in  the  writings  of  R.  C.  Jones.  This  is  to  use  the  form 


W«(k  ,  k  ) 
B  X  Y 


M 


a 


I  Z 
k  +  k 
X  y  _ 


1-3^ 


(21) 


If  ct  is  not  too  large,  this  preserves  agreement  with  the  experi¬ 
mental  results  as  out  subsequent  discussion  will  show. 

First,  let  us  obtain  the  correlation  function  associated  with  Equa¬ 
tion  (21);  we  must  evaluate  the  Fourier  integral, 


/oo  ^oo  r  2 

[  M  lo  +  k 

CO  CO  L 


k 

X 


1 


-3/2 


y  j 


exp  k^x  +  k^yTJ  dk^  dk^ 


(22) 


i  1 


This  turn*  out  to  be 


where 


fix.  y) 


itM  ^-2iTod 
a 


d 


2  ^ 

+  y2) 


(23) 
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This  result  shows  that^tx,  y)  is  a  function  only  of  d=(x^  +  y^)*' 
and  thus  indicates  that  the  choice  of  the  form  (21)  is  consistant  with  our 
previous  assumption  of  isotropic  background.  This  property  of  the  cor* 
relation  function  (jdependence  upon  (x^  +  y^)^3  can  be  proven  to  hold  for 
power  spectra  of  the  form 


W(k^. 


function  of  (k 


2 

X 


y 


(24) 


The  measurements  of  W'(k^,  k^)  that  are  known  to  the  writer  really 
measure  the  one- dimensional  power  spectrum  —  more  precisely,  the  one- 
dimensional  Fourier  transform  of  die  correlation  function.  For  the  form 
(23)»  this  one- dimensional  transform  is 


W  ,k,  =  r"  ^  ^-Z-idk 

®  1“! 

M  1 

4v|a(  ^2  + 


(25) 


2  2,  ^ 

The  function  (a  h  )  ,  a  normalized  form  of  Equation  (25),  is 

plotted  in  Figure  2,  versus  k  for  various  values  of  a-  It  is  seen  that 
asymptotically  for  high  wave  numbers,  this  spectrum  follows  a  k  law. 
The  value  of  a  used  determines  how  the  curve  deviates  from  the  (-2) 
law.  The  several  sets  of  measurements  results  that  have  been  seen  by 
ttie  writer  (schematically  indicated  by  dashed  line  in  Figure  2),  are  not 
iacoosistent  with  this  type  of  function,  but  it  is  difficult  to  make  a 
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I  iO  100  1000 


k  (WAVES  PER  RADIAN) 

FIG.  2  THEORETICAL  ONE- DIMENSIONAL  POWER  SPECTRUM. 
DASHED  CURVE  INOIC^ATES  APPROXIMATE  SHAPE  OF 
MEASURED  POWER  SPECTRA  (LEVEL  ARB'TRARY). 


jud^fement  ae  to  what  the  best  value  of  a  would  be.  Perhaps  various 
background  conditions  (which  give  widely  differing  absolute  levels  of 
background)  all  give  the  same  asymptotic  {  Z)  law,  but  vary  in  their 
effective  values  of  a>  If  all  that  is  known  about  the  measured  spectrum 
iSi  for  esample^  that  it  follows  a  fairly  accurate  (*2)  law  in  the  range 
say  above  10  waves  per  radian,  then  it  is  seen  that  any  velue  of  a  up  to 
about  S  is  consistent  with  such  measurements.  This  suggests  that  since 
the  high  k  part  of  the  spectrum  has  a  form  independent  of  some  of 
the  results  of  the  subsequent  analysis  of  output  noise  may  also  be  inde¬ 
pendent  of  cr»  As  measurements  become  more  complete,  it  may  be  pos¬ 
sible  to  choose  more  accurately  between  values  of  a  or  perhaps,  to  pick 
a  more  accurate  form  for  the  autocorrelation  function.  It  is  not  unlikely 
tiiat  as  measurements  of  background  become  more  accurate,  the  (-Z)  law 
and  the  form  e  may  have  to  be  modified*  Measurements  aimed 

at  directly  measuring  the  background  correlation  function  could  perhaps 
settle  this  point  more  quickly  than  the  spectral  measurements  For  ex¬ 
ploratory  purposes,  however,  it  seems  to  be  of  some  value  to  study 
Gin)  and  its  Fourier  transform  for  a  range  of  values  of  a-  At  the  time  of 

writing*  the  function  G(u)  and  its  transform  are  being  numerically  studied 

2 

For  the  time  being,  the  Gaussian  form  €xp{-pd  )  for  the  correla¬ 
tion  function,  which  is  often  a  likely  candidate  in  similar  problems  has 
been  eliminated  from  consideration  because  the  range  of  the  effective  (-d) 
law,  which  seems  to  be  an  octave  or  more  in  the  known  measurements, 
cannot  be  fitted  very  well  to  the  transform  of  exp(-pd^)  which  is  of  the 
form  expi-qk^);  p  and  q  are  constants  in  the  above;  d  and  k  have  signi¬ 
ficance  similar  to  that  in  Equations  (22)  and  (25) 


STUDY  OF  A  PARTICULAR  BACKGROUND  CORRELATION  FUNCTION 

Stu<lies  and  computations  are  under  way  of  the  properties  of  G(u) 
and  its  transform  W^(z).  At  the  time  of  writing  these  studies  are  in- 
Caesq^late,  but  a  few  preliminary  remarks  may  be  made. 

Substituting  the  form  (23)  into  the  expression  for  G(u),  we  have 


C5<u) 


^  r  r\r' 

®  Jo  Jo 


(26) 


1  3 


dr  dr 


Introducing  the  dimensionless  parameters 


r 

a 


ZttaSL  ~  h  , 


{27) 


we  obtain 


Fh{«) 


G(«) 


^  a 


r  r 

Vo  Vo 


dx^Xj 


exp 


h  (x*  +  X  “  -  2x 


2 

2 


J  x^  cos  u) 


i/Z 


(28) 


where 


C 


IM 

a 


Fj^(u)  is  the  significant  function  to  be  studied.  The  parameter  h  may 
be  written  as 


■  T/z^  • 

Now  i/Zita  may  be  interpreted  as  a  "correlation  length,  "  that  is,  l/Zva  is 
the  value  of  distance  d  for  which  the  (normalized)  correlation  function 
e  falls  to  one  e^  of  its  peak  value  (which  occurs  for  d=  0).  The 

greater  l/2xo,  the  greater  is  the  relative  correlation  of  two  points  sep¬ 
arated  by  a  distance  d.  That  is,  the  normalized  correlation  function  falls 
more  slowly  from  its  peak  value  for  high  values  of  the  correlation  length 
1/2x0-.  The  parameter  h  is  thus  the  ratio  of  the  radius  of  the  field  of 
view  and  the  correlation  length  of  the  background.  The  significant  quan¬ 
tity  affecting  the  properties  of  F^{u)  is  seen  from  Equation  (28)  to  be 

this  dimensionless  parameter  h.  In  addition,  there  is  a  scale  factor  of 
4 

a  ,  which  simply  means  that  all  other  things  being  equal,  the  output  noise 
has  a  power  level  which  is  proportional  to  the  square  of  the  area  of  the 
field  of  view. 

The  dependence  of  F^(u)  upon  h  and  u  is  sketched  in  a  purely 
qualitative  way  in  Figure  3.  The  range  0<u^2x  is  shown;  F^(u)  is 
periodic  in  u  (a  result  of  the  assumption  v=  0)  The  behavior  of  the 
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transform  of  very  roughly  as  follows:  The  higher  values  of 

h  will  yield  a  slower  high  frequency  faii-off  (more  relative  high  frequency 
content);  but  the  range  and  extent  of  these  variations  can  only  be  obtained 
from  the  complete  study. 

To  obtain  a  rough  idea  of  what  values  of  h  may  be  of  interest,  let  us 
imagine  that  the  condition  referred  to  above  holds,  namely  that  the  meas¬ 
ured  spectrum  follows  roughly  a  (-2)  power  law  in  the  range  above  tO 
waves  per  radian;  in  the  absence  of  further  experimental  details,  any  value 
of  o  up  to  about  0=5,  may  provide  a  satisfactory  fit.  JLet  us  take  the 
radius  a  to  be  2  degrees.  Then  h  =  2ttaa  =  2w5(2tt/180)*  l.i 
The  numerical  work  will  cover  values  of  h  from  0  to  5. 

Once  1*  computed  for  a  sufficient  range  of  h,  it  will  serv^e  the 

purpose  of  a  '’universal  curve"  applicable  to  any  value  of  a  and  any  reticle 
function  R(e).  Its  transform  will  likewise  represent  the  power  spectrum 
of  the  noise  as  seen  by  any  reticle  function.  By  multiplying  this  transform 
by  the  square  of  the  transform  of  the  reticle  function,  Equation  (16),  we 
nriay  obtain  the  spectrum  of  the  noise  output  of  the  reticle.  This  may  then 
easily  be  combined  with  transforms  of  following  stages. 

This  numerical  work  is  only  valid  for  the  background  correlation 
function  e  jhe  method  of  attack,  however,  is  still  valid  for  any 

improved  background  functions  that  may  later  be  obtained. 
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APPENDIX  A 


TRANSFORMATION  OF  A 
CERTAIN  INTEGRAL  EXPRi^lON 


This  appendix  concerns  itself  with  the  transformation  of 


Jo  Jo 


AB'  G(©'  -  6)  R<0)  R(8*  -  ut)  . 


wiiexe  G(x)  s  G(x  *  Zn) 
and  R(x)  =  R(x  -f  Zv) 

Xatxodttce  the  change  of  varialdes 
u  a  0'  -  «T  -  0  , 


V  a  6'  <-  ta»T  -f  0 


■  _  VI  +  V  , 

0  -  ■"  — — '  -f  U)T 


-  V  -  u  , 

0  =  . .  4  U>T 


The  Jacobian  of  this  transformation  is 


1/2  -  1/2 


1/2  1/21 


=  1/2 


We  tiien  have  (see  Figure  A-1  for  geometrical  interpretation) 


21=  f  duG(u  +  oit)  f  dv  ^ 

J-Z^  -WT  •' V  =  2wt  -  u  \  Z  I  \  Z  j 

/•  2v  -  <*>T  r  -ZuiT  -  u  /  \  /  . 

+  j  duG(u  +  <*»t)  J  dv  r(  — - ~]^(~ - - 

J‘0»T  u  \  2  /  \  2 


letroduciag  the  change  of  variable  ^  ^  =  s  we  have 


FIG.  A-1 


geometrical  interpretation  of  CHANi 
OF  VARIABLE  IN  TEXT 


2  jdB  R{8)  R(B  +  U)  . 


•o  that 


■/, 


-«T  r Zv 

1  *  i  «iu  G(u  +  wt)  f  ds  R{s)  Rfs  +  u) 

•  Ztr~ur  ■ 


■L 


WT-U 

2w-u>t  r  iw-wT-u 

du  G(u  +  «t)  j  ds  R(8)  R(s  +  u) 

«T  ^  0 


Now 


iT  2* 

-/-WT-a 


2*  ^2it 

d*  R{8)  R(b  +  u)  =  /  ds  R(s)  R(s  +  u  +  2tr) 

-u)t-u-2it 


^  -tiJT“U 
✓  -WT-u-Zir 


ds  R(a)  R(s  4  u  -f  2ir)  . 


Calling 


^  2w 

I  ds  R{8}  R(a  +  u)  ds  =  P(u)  =  P(u  +  2ir)  , 

y  0 


we  have 


/-«T  r  ^  Zir 

du  C<u  +  wt)  j 

2Tr-«T  -WT-U 


ds  R(8)  R(s  +  u  +  2ir)-P(u  +  2ir)] 


■2tT 


.  r 


-«T 


■/_ 

■I. 


^  2ir-u>T 
-2Tr-uT 


du  G(u  +  2tr  +  ujt) 


dv  G(v  +  u)t) 


wr 

-U)T+  2ir 
wr 


WT-U-j 

/Ztt 

ds  R( 

-  -WT-V 


ds  R(8)R(8  +  u  +  2Tr)-P(u  +  Ztt) 


8)R(s  +  v)-P(v) 


du  G(u  +  wt) 


f 


WT 

Zrr 


R(s)R(s  +  u)  ds  -  P(u) 


-WT 
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so  that 


where 


r  2s 

I  =  /  du  G(u  4  wt) 


•«T 


ir” 

Lr  "WT-u 


•)RC»  -i-  u)  ds 


2ir-«T-u 


R{a)R(8  +  u)da  -  P«u) 


_  r 


■«T+  2ir 


[/:: 


=  I  du  G(u  +  u>t)  1  j  R(8)R(8  ■¥  u)d«  -  P{u) 

J~u»T  L-'-Zir 


-WT 


-«T+  Zit 


du  G<u  =  «t)  P{u) 


P(u)  =  /  R{8)  R(p  +  u)  da  . 

•^0 
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APPENDIX  B 


ANALYSIS  OF  BACKGROUND  MEASUREMENTS 
WITH  CIRCULAR  SCAN 


In  coKonection  with  the  study  of  the  effect  of  rotating  reticle  systems 
background,  th.ere  is  at  present  an  interest  in  a  program  of  back¬ 
ground  '*gracltent"  measurements  in  which  the  background  is  scanned  with 
a  smaU  field  of  view  moving  repetitively  around  a  circle  of  radius  a.  The 
Fooriei  coefficients  of  the  periodic  output  are  taken  with  the  object  of 
relating  these  Fourier  coefficients  to  the  properties  of  the  background. 
C^orating  under  the  assumption  that  the  background  is  a  stationary  ran- 
dmm  fiinction,  we  give  an  indication  of  the  kind  of  reductions  needed  to 
convert  such  data  into  a  background  power  spectrum  or  correlation  f^c- 
lioa. 

EXPRESSIONS  FOR  GF,NERAL  PATHS 

Before  specializing  to  the  case  of  the  circular  path  we  present 
acsne  expressions  which  are  valid  for  general  paths. 

Let  us  assume  that  our  instrument  measures  the  radiance  B{s) 

M  a  fimction  of  the  running  angular  coordinate  s  along  the  path,  and 
fitis  type  of  scan  is  repeated  over  many  similar  segments  of  the 
field,  and  that  for  each  measurement  along  the  path  we  obtain  a  set  of 
Jfottrier  coefficients. 

A  typical  complex  Fourier  coefficient  is 

.  -i2ir(n/A)s  , 

B(8}  e  '  '  '  ds 

fdlMre  A  is  the  length  of  the  path. 


footnote  (1)  of  the  main  text. 


(1) 
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The  square  of  the  magnitude  of  C{n)  is 


C{n)  C*(n)  = 


^-i2Tr(n/AHs-p)  ^ 


(2) 


where  n  is  any  integer  between  >a>  and  oo . 

We  shall  assume  that  the  measurements  give  us  tfee  C^n)  C*^a)  . 
To  obtain  average  data  we  average  over  many  paths  or  cases. 
We  have 


=  <^C(n)  C*(n) 


^  ^-i2TT(n/A)(8-p)  ^ 


j  J 

«'0  -'O 


B{s)  B(p)^  ds  dp 


(3) 


y  ^ 

where^  ^denotes  the  average  over  many  cases  (stochastic  or  ensem- 


ble  average). 


RESTRICTION  TO  ISOTROPIC  BACKGROUND  AND 
CIRCULAR  OR  STRAIGHT- UNE  PATHS 


Further  evaluation  of  Equation  (3)  depends  on  the  assuunptions 
that  are  made  on  the  form  of  <^^^(6)  B{p)^  .  Let  us  assume  for  the 
moment  that 


B(8)  B(p) 


p(8  -  p)  “  P(p  -  s) 


<4) 


i.  e.  ,  that  the  average  value  of  B(8)  B(p)  depends  only  upon  the 
difference  between  the  rimning  coordinates  s  and  p.  It  will  be  seen 
below  that  this  assumption,  which  enables  us  to  make  progress  iu  the 
interpretation  of  Equation  (3),  implies  certain  restrictions  on  the  back¬ 
ground  and  on  the  paths. 
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Withthi*  astumptioa.  Equation  (3)  becomes 


Jq  Ja 


D  = 

®  a‘  •'0  -'o 


Xf  now  change  to  the  variables 


/A)(b-p) 


^{s  -  p)  ds  dp 


(5) 


U  =  8  -  p 
V  =  8  +  p  , 


w«  Obtain,  noting  that  the  Jacobian  of  this  transformation  is  l/2 , 


D  a 
n 


r  dn  r  ^  0(u) 

a2  <A  “'IuI 

^  (2A  -  2jul)  ^  du 


JL 

2a2  ‘'-a 


A  -A  \  A  / 


C08  ^(u)  du 


Thus  with  the  help  of  the  assumption  (4)  we  have  reduced  (3)  from 
a  doable  integral  to  a  single  integral.  The  interpretation  of  this  inte- 
will  be  discussed  below;  in  particular  we  shall  examine  the  extent 
|0  which  information  about  the  background  random  hmction  can  be  ex- 
tvactad  from  (8).  fTrst  however,  we  should  pause  to  consider  under 
Oddch  ccmditions  ^B(8)  B(p)  y  will  be  of  the  form  ^(s  -  p),  so  that  the 
fnrm  (8)  will  be  valid.  Let  us  assume  that  the  two-dimensional  noise 
fa  isotropic;  this  means,  as  is  discussed  above,  that  the  correlation 
JtoCtion»  which  is  with  full  generality  a  function  of  two  space  variables. 


(6) 


(7) 


(8) 
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say  the  orthogonal  variables  x  and  y,  is  in  this  case  a  function  only  of 
d  =  [^x  +  y  J  In  other  words  <^B(8)B(p]^,  which  is  a  correlation 

function,  is  a  function  f(d)  of  the  scalar  distance  between  the  path 
points  labelled  s  and  p.  However,  this  is  still  not  sufiicient  for  the 
validity  of  {4),  in  which  ^B{s)Blp^  is  set  equal  to  a  function  of  s  ~  p. 

In  order  for  this  to  be  valid,  the  distance  between  the  points  s  and  p 
should  depend  only  upon  the  quantity  s  -  p,  and  not  upon  s;  this  condi> 
tion  seems  to  be  satisfied  only  if  the  scanning  path  is  a  straight  line  or 
a  part  of  a  circle.  For  an  ellipse  for  example,  the  same  path  length 
difference  (s  -  p)  will  give  different  distances  between  s  and  p,  depend¬ 
ing  upon  what  part  of  the  elliptical  path  s  and  p  refer  to. 

Thus,  formula  (8)  refers  to  the  case  of  isotropic  noise  and  straight 
line  or  circular  scans. 

Equation  (8)  exhibits  as  the  cosine  transform  of  the  quantity 
{1  -  u/A)  ^{u).  A  possible  procedure  then,  is  to  invert  (8)  to  obtain 
(1  -  u/A)  ^(u),  divide  by  1  -  u/A  to  find  ^(u)  and  then  from  ^(u)  find 
)^(d),  the  correlation  function  of  the  backgroimd.  If  the  path  s  has  been 
a  segment  of  a  straight  line,  then  u  is  the  same  as  d  and  ^(u)  =  ^(d)=  V* (d). 
If  the  path  s  is  a  circular  path  (or  part  of  a  circular  path),  then  the 
relation  between  u  and  d  is  that  d  is  the  chord  connecting  the  ends  of 
the  circular  arc  whose  length  is  u.  That  is,  if  the  radius  of  the  circle 
is  a,  (d/Za)  =  sin  ^u/Za).  We  have  then 


'y^(d)  =  0[Za  sin  (d/Za)]  (9) 

so  that  by  a  straightforward  nonlinear  scale  change  we  can  obtain  Y'  (d) 
having  once  obtained  ^(u).  We  may  then  use  Y'  (d)  to  obtain  G(u)  the 
noise  correlation  function  by  means  of  Equation  (9a)  of  the  main  body 
of  this  paper. 

The  main  task  then,  is  the  inversion  of  Equation  (8)  to  obtain 
(1  •  u/A)  p(u).  We  have  noted  that  the  is  expressed  as  a  cosine  trans¬ 
form  of  (1  *  u/A)  4(u).  More  precisely,  however,  the  cosine  transform 
of  (1  -  j  u|  /A)  9(u)  that  fully  represents  the  function  in  the  range  o  to  A  is 
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E  » 
n 


0(u)  cos 


2iTnu 


ixiag  Equation  (iO)  with  Equation  (8)  we  see  that  the  D  are  not 

i4«ittical  with  the  £  but  that  E  ^  D  /^.  Since  however  the  D 

n  n  n/  2  n 

vailable  only  for  integral  n,  half  the  £^  are  missing.  To  obtain 
£  we  must  interpolate  between  the  values  of  D  ;  that  is  for 
le  on  tihe  curve  of  versus  n,  E^  is 

Sauce  having  obtained  a  complete  set  of  E^,  half  of  which  are  the 
i  the  others  interpolated  values,  we  may  find  (1  -  u/A)  ^(u)  by  the 
Um  formula 


(-i) 


^(u)  = 


£  cos 
n 


rile  task  of  interpolation  between  the  may  be  rendered  some- 
lore  hopeful  by  the  fact  that  (1  -  u/A)^(u)  should  be  a  correlation 
m  eo  that  none  of  the  £^  should  be  negative.  Having  obtained 
'A|#<u}  the  rest  of  the  procedure  outlined  above  for  obtaining  ^ (d) 

I  followed.  If  desired  one  may  then  obtain  the  one -dimensional 
spectrum  of  the  background  by  Fourier-transforming  S^(d). 

It  ie  seen  from  the  above  discussion  that  the  procedure  of  obtain¬ 
ed)  from  the  D  would  be  subiect  to  considerable  error.  This 
n 

tt3  if  our  aim  lo  :ompute  G(u),  iu  which  /  (d)  itself  appears, 
m  measurements  of  background  be  arranged  to  directly  obtain  (d) 

than  data  like  the  O  . 

n 

ritue  for  example,  one  could  measure  the  average  value  over  s 
)  |l(e  ^  u)  where  s  and  s  +  u  are,  say,  positions  in  the  circular 
eth*  This  average  value,  itself  averaged  over  a  sufficiently 
Mmber  of  cases,  is  then  the  function  ^(u)  discussed  above. 
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y^id)  is  then  easil)*^  obtained.  Straight  line  paths  could,  of  course, 
also  be  used,  yielding  y^(d)  directly. 

A  further  remark  needs  to  be  made  concerning  the  D^.  It  can  be 
shown  ^  that  the  integral 

\  p(u)  cos  2Trfu  du 

A  / 

which  appears  in  Equation  (8j|,  approaches  Ztrfu  du  as  A 

approaches  infinity.  This  last  integral,  however,  is  die  power  spectrum 
associated  w’.th  the  correlation  fiinction  p(u).  As  A  approaches  infinity, 
moreover,  p{u)  approaches  y^{d).  Therefore  for  large  A  the  are 
approximately  equal  to  ^l/A)  {n/A),  where  the  one  dimensional  power 

spectrum  (kj  is  defined,  as  in  Equation  (25)  of  the  main  text,  as  the 
Fourier  transform  of  ^(d): 


Wg  (k) 


V" 


(d)  e 


-2iridk 


dd 


yoo 

=  /  "y/"  (d)  cos  2-Trdk  dd  . 


For  a  given  value  of  A,  the  statement  is  more  nearly  true,  on  the  aver¬ 
age,  for  high  n. 


Titchrr.arsh. 


"Theory  of  Fourier  Integrals." 


p.  36. 
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